Confronting General Relativity with Further Cosmological Data 
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Deviations from general relativity in order to explain cosmic acceleration generically have both 
time and scale dependent signatures in cosmological data. We extend our previous work by investi- 
gating model independent gravitational deviations in bins of redshift and length scale, by incorpo- 
rating further cosmological probes such as temperature-galaxy and galaxy-galaxy cross-correlations, 
and by examining correlations between deviations. Markov Chain Monte Carlo likelihood analysis 
of the model independent parameters fitting current data indicates that at low redshift general rela- 
tivity deviates from the best fit at the 99% confidence level. We trace this to two different properties 
of the CFHTLS weak lensing data set and demonstrate that COSMOS weak lensing data does not 
show such deviation. Upcoming galaxy survey data will greatly improve the ability to test time and 
scale dependent extensions to gravity and we calculate the constraints that the BigBOSS galaxy 
redshift survey could enable. 



I. INTRODUCTION 

Gravitation is the key force governing the expansion 
and evolution of the universe. The unexpected obser- 
vations of cosmic acceleration may indicate that some 
aspects of this fundamental force remain a mystery. Gen- 
eral relativity is a hugely successful theory of gravity over 
the ranges it has been tested, but we should continue to 
test it in greater detail in regions, such as on cosmic 
scales, where it has not been sufficiently probed. 

Since it is not clear what form deviations from general 
relativity (GR) may take, it is useful not only to adopt 
specific models extending GR but also to consider model 
independent approaches. These generally parameterize 
the relation between the metric potentials, the relation 
between the potentials and the matter density, or similar 
forms. A translation table between many of the most 
common conventions was provided in 

The effects giving rise to cosmic acceleration must take 
place on the largest length scales, but general relativity 
is known to be highly accurate on small scales (solar sys- 
tem and laboratory), so the deviations must have scale 
dependence. This can either be innate (from the scale de- 
pendence in the Poisson equation), or explicit. Similarly, 
conditions in the early universe such as during primordial 
nucleosynthesis or recombination can be well explained 
within GR, and acceleration is a recent phenomenon, so 
the deviation from GR should also be time dependent. 

In this article we broaden consideration of the devia- 
tion parametrization, and convert to more observation- 
ally direct variables than used in [l|. In Sec. [TT] we ex- 
amine some possible time and space dependencies and 
examine the correlation between the deviation variables. 
By adopting a model independent, binned formalism we 
avoid putting in ad hoc assumptions about the form of 
the deviation, letting the data determine the results. We 
consider different data types probing the matter den- 
sity distribution in Sec. IIII1 going beyond the cosmic 
microwave background (CMB) perturbations, Type la 



supernova distance-redshift relation, and weak gravita- 
tional lensing used in [lj. Prospects for further improve- 
ments in constraints from future data are investigated in 
Sec. EH 



CONSTRAINING DEVIATIONS OF 
GRAVITY 



II. 



The relations between the two metric potentials (often 
called the gravitational slip), the matter density and ve- 
locity fields (continuity equation), the matter density and 
a metric potential (Poisson equation), and velocity field 
and the other metric potential (Euler equation) form a 
system of equations describing the spacetime and its con- 
tents. Modifications to gravity adjust these interrelations 
and so one can parameterize these theories by inserting 
time and space dependent functions in the usual GR re- 
lations. 

One example is to define the gravitational slip as 



ip = [1 + w(a, k)] <f> . 



(1) 



where the metric is given in conformal Newtonian gauge 
through 



ds 2 = a 2 [-(l + 2ip) dr 2 + (1 - 2<f>) dx 2 ) 



(2) 



and a is the scale factor, k the wavenumber, r the con- 
formal time, and x the spatial coordinate. Preserving 
stress energy conservation, and so the continuity and Eu- 
ler equations, we are left with needing the Poisson equa- 
tion, modified to 



— k 2 <fi = ^ttGno 2 p m A-m x M a > ■ 



(3) 



where Gn is Newton's constant, p m is the homogeneous 
part of the matter density and A m the perturbed part 
written in gauge-invariant form, i.e. 



ZHO r , 
k 2 



(4) 



2 



in the notation of , where 9 m is the velocity perturba- 
tion and % is the conformal Hubble parameter. 

These two functions, w and fx, were used in [l[ and 
several other papers, and are equivalent to many other 
parametrizations as detailed in the translation table of 
In this paper, we will present a few further results 
using these variables, but the bulk of the paper will use 
"decorrelated" parameters based on these. 



A. zu and fi as Variables 

One of the main focuses in [l[ was to test consistency 
with GR. For this, only one of w or fi were varied at 
a time. Since a shift in w could be compensated by a 
corresponding shift in \x (see Fig. 2 of [l], or the degen- 
eracy line in their Fig. 7) to preserve the observational 
agreement, the intent of varying one at a time was to 
make it more difficult to achieve agreement with GR and 
hence provide a more conservative result. Despite this 
"handicapping" , agreement with GR indeed occurred. 

In Fig. [T]we show what happens when fi is allowed to 
vary simultaneously with vu, treating both functions as 
being composed of constant values within each of three 
rcdshift bins (z < 1, 1 < z < 2, and 2 < z < 9, with z > 9 
fixed to GR) . To compare tollf , we use data constraints 
from WMAP 5 year CMB Union2 supernova dis- 

tances [f|, and COSMOS weak lensing 0] data sets. The 
narrow ID distribution of w recreates the fixed /x case of 
the middle panel of Fig. 5 of while the wider distri- 
bution shows the results for w when also fitting \i. The 
68% cl range increases by approximately a factor of 6. 
The other rcdshift bins behave similarly. Thus, the accu- 
racy of measurement of the deviations w (and /z) is not 
particularly tight. 

The degeneracy between the two post-GR functions is 
clearly seen in the 2D probability distributions of Fig. [2] 
The banana shape discussed in [l| persists here, even 
though we use independent bins of redshift rather than 
the a functional dependence assumed in their Fig. 7. 
The solid black curve shows a theoretically motivated 
compensation relation largely responsible for the degen- 
eracy. 

We will be able to improve the constraints, and our 
understanding of the observational leverage on modified 
gravity, by "trading" precision on one combination of w 
and n for that of another combination. This basically 
corresponds to choosing variables along and perpendicu- 
lar to the main degeneracy direction, as we now discuss. 



B. Separating Parameter Effects: Q and V 

Several types of observables are predominantly sen- 
sitive to the sum of the potentials, e.g. the integrated 
Sachs- Wolfe (ISW) effect and gravitational lensing. Writ- 
ing the Poisson equation in terms of such a sum yields 




vary only C3 
vary both C5 and |l 



G3 (l<z<2) 



FIG. 1: ID marginalized probability of the post-GR param- 
eter zu in the redshift bin 1 < z < 2. The narrower, dashed 
(black) distribution fixes the other post-GR parameter /i = f , 
making it more difficult to fit the data (consistent with GR) 
by compensating one parameter with another. The wider, 
solid (red) distribution includes a simultaneous fit for fi. All 
other cosmological parameters, including zu and [i in the lower 
and higher redshift bins, are marginalized over. 



=i- 




k 2 ((j) + ip) = 8irG N a 2 p m A m x Q , 



(5) 



FIG. 2: 2D joint probability contours at 95% cl between the 
post-GR functions zu and n, for three independent redshift 
bins. Values within the redshift bins are consistent with each 
other and with GR (denoted by the cross at (0,0)). The solid, 
black curve, motivated by a modified Poisson equation, closely 
follows the degeneracy direction, and suggests a more insight- 
ful parametrization using variables along and perpendicular 
to the curve. 



where Eqs. (p} and show that Q = [i (2 + zu)/2. It is 
not surprising therefore that the confidence contours in 
the [i-m plane are banana shaped with strong curvature. 

The main degeneracy curve illustrated in Fig. [5] is pre- 
cisely the combination entering Q. Therefore it makes 
sense to switch variables to use this combination as one 
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parameter. We can also use Eqs. §5§ and ((3j) to define a 
Poisson-like equation for tf>; here we write all three equa- 
tions together to show the parallelism: 

- fc 2 (0 + V) = 8nG N a 2 p m A m x Q (6) 
-k 2 4> = AirG N a 2 p m A m x /i (7) 
~k 2 %j) = 4:TTG N a 2 p m A m x V . (8) 

The parameter V is precisely the parameter identified as 
mostly sensitive to growth of structure in [l[ (there called 
S; note that [8| earlier noted this set to be of interest, 
calling V as fj,, and 5 as S; @ also explored this later). 
It is also closely related to the growth index parameter 7 

nana. 

The new post-GR functions are related to the old ones 
via 



Q = li- 

V = + ro) 



At = 2£ - V 

2V- 20 
2Q-V 



(9) 
(10) 



We will see that the new functions are substantially 
decorrelated from each other, producing more indepen- 
dent constraints when using the observational data. (The 
symbol Q is meant to evoke an effective Newton's con- 
stant in the total Poisson equation; V recalls the "veloc- 
ity" equation arising from the relation between the po- 
tential rp and the matter velocity field, central to growth 
of structure.) 

One expects that the integrated Sachs- Wolfe effect 
and, in large part, weak gravitational lensing data will 
mostly constrain Q (i.e. across the degeneracy direction 
seen in Fig. [2]) and have little leverage on V (i.e. along 
the degeneracy direction seen in Fig. [5]). Probes that 
involve growth, such as weak gravitational lensing and 
the cross-correlation between the CMB and the galaxy 
density field, to some extent, and the galaxy-galaxy den- 
sity power spectrum, should place some constraint on V. 
From Fig. [5] one expects that current weak lensing data 
will not be that strong, however, so we will also inves- 
tigate the role of current density field data in Sec. IIIII 
Future galaxy survey data should tighten the constraints 
further; see Sec. HVIfor further discussion. 



C. Redshift and Scale Dependence 

The post-GR functions will generally be functions of 
both time (redshift) and length scale (wavenumber) . Wc 
do not necessarily want to assume a particular functional 
form, so we begin by allowing the values of Q and V to 
take arbitrary values within independent bins of redshift 
z and wavemode k (also see early work by [Hjl)- 

If we examine the redshift dependence of G, using no 
scale dependence initially, we find that the values of Q in 
different bins are positively correlated. We consider two 
independent redshift bins, with < z < 1 and 1 < z < 2. 
For z > 2 we assume the GR values. The characteristics 



discussed below do not change if we add a third bin at 
2 < z < 9, but the remainder of this paper uses two bins. 
The degeneracy direction between Q(0 < z < 1) and 
Q(l < z < 2) corresponds roughly to a dependence on 



scale factor G(a) 



at least for z < 2. The function V 



shows a negative correlation between V(0 < z < 1) and 
V(l < z < 2), such that they roughly compensate each 
other: V(0<Kl)-l« -[V(l < z < 2) - 1]. 

Regarding degeneracies with other cosmological pa- 
rameters, there is little correlation except with the mass 
fluctuation amplitude as . This accords with the principal 
influence of V and Q being on growth of scalar pertur- 
bations, especially at late times. The main effect is a 
positive correlation between V(l < z < 2) and as; recall 
that V is the post-GR parameter most strongly entering 
into the growth of A m . That the higher z bin of V is most 
correlated follows from growth being cumulative, so the 
higher redshift bin has a longer lever arm of influence 
to imprint the effects of gravitational modifications. We 
also find a slight negative correlation between Q and as- 
This is related to the weak lensing data, which involves 
the sum of the potentials as well as the growth (see the 
discussion at the end of Section |TTT]) . For higher Q, lower 
values of as will produce comparable lensing potentials. 
Thus, larger Q does not cause as to decrease per se (the 
way larger V amplifies growth), rather it brings lower 
values of as into better agreement with the data. 

Now considering scale dependence, we introduce two 
bins in wavenumber k, running from k = 10~ 4 — 
1CT 2 Mpc -1 and k > 1CT 2 Mpc -1 . The low k range 
represents the large scales from roughly Hubble scale to 
matter-radiation equality horizon scale, and the high k 
range corresponds to scales roughly over which non-CMB 
probes have leverage. For example we expect that the 
matter power spectrum (including weak lensing) would 
mostly constrain the second bin. Thus in total we fit for 
8 post-GR parameters: Q and V values, each in 2 bins of 
z and 2 bins of k. 

In Fig. [3] we plot the 68% and 95% confidence limit 
contours in {V, G} space for all bins of z and k. These 
contours have been calculated generalizing the modified 
COSMOMC code used in ^ and incorporate WMAP7 
13 1, supernova Union2 Q, and CFHTLS weak lensing 
14 1 data. The original COSMOMC was presented in 
and the weak lensing data module is from [18j. 

From Figs. [3(a)p(c) we see that our initial supposition 
that V and Q are mostly independent (or, at least, less 
correlated than zu and n are in Fig. [2]) is correct. We also 
see that the constraint on Q is in all cases stronger than 
the constraint on V. 

The dot-dashed line on the right side of the figures 
corresponds to values of V and Q for which [x = (see 
Eq. [3]). This would imply that the metric is independent 
of matter perturbations, which seems unphysical. For the 
most part, restriction to fi > docs not strongly affect 
the contours. However, the data considered so far is not 
so strong as to exclude the fi < region (to the right of 
the line) without imposing a prior. We feel the prior is 
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FIG. 3: 68% and 95% confidence limit contours for V — 1 and Q — 1 are given for 2x2 binning in redshift and k space, using 
WMAP7 Union2 Q, and CFHTLS Q data. The diagonal, dot-dashed line denotes values of V and Q for which fi = 
and gravity vanishes (see Eq. [3}. The x's denote GR values. 



justified in that, referring back to Eq. (A5) of \i < 
implies (as a consequence of stress-energy conservation) 
that there will be some value of k for which a factor on 
the left hand side of that equation goes to zero, causing 
A m to diverge. 

The upper two figures show the results for the low k 
bin, where the current data is most constraining. Note 
that GR is comfortably within the 68% cl contour. The 
constraints above z « 1 arc slightly tighter, since the 
ISW effect is more sensitive to this region. The bottom 
two figures give the results for the high k bin, and here Q 
is significantly more constrained at higher z, again due to 
the ISW. Comparison of the different k bins in the same 
redshift range show that low k is better constrained, due 
to the ISW. 

Figure |3(d)| the high k - low z case, exhibits a num- 
ber of peculiarities. It is much less constrained than the 



other cases, and shows a higher correlation between Q 
and V. There is also an apparent, nearly 3cr exclusion 
of General Relativity. Although this is a tantalizing re- 
sult, it should not be taken too seriously. Since the ISW 
effect is an integral over redshift, the low redshift bins 
have very weak effects on the CMB anisotropy spectra 
and cannot be tightly constrained by WMAP7. There- 
fore, any systematic errors that create tension between 
the CFHTLS data and WMAP will be able to manifest 
themselves as non-GR values of V and Q in these bins. 
Furthermore, the high k bins encompass scales for which 
the ISW effect is subdominant anyway. It behooves us to 
turn our attention, then, to other data sets that may be 
sensitive to these bins in the hope of strengthening our 
confidence in these constraints. We do this in the next 
section. 

Before we proceed, however, we note that this is not 
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the first work to find a 2a exclusion of GR at small scales 
and low redshift. Reference @ reports a similar result 
for their parametrization \n ( see Section IVB of that 
work). One curious difference, though, is that their anal- 
ysis shows a preference for Q < 1 (their E is equivalent 
to our Q), whereas our Fig. |3(d)"| shows a clear preference 
for Q > 1. This difference can be traced to the different 
k binning schemes. The discussion in Q attributes the 
preference for non-GR Q as a means to fit a systematic 
bump in the CFHTLS weak lensing data at large scales 
(see their Fig. 8; they state that the CFHTLS team as- 
cribes this to residual systematics; also see Sec. 4.3 of 
flip. They divide their k bins at k = 0.1 h Mpc -1 . This 
is approximately where the bump in the CFHTLS data 
occurs (for z < 1). The MCMC code exploits this by 
selecting a large value of TO to increase the overall lens- 
ing amplitude and fit the bump at large scales (and low 
k) while reducing the value of Q in the small scale (large 
k) bin to prevent that increased amplitude from spoiling 
the fit to the smaller scale data. This allows them to 
alter the shape of the weak lensing power spectrum to 
rise and fall with the data. Since we divide the k bins at 
k = 0.01 Mpc -1 , however, the same shift in parameters 
would suppress growth, and hence weak lensing power, 
over too large a range of angles. 

Figure E] illustrates this, as well as recreating Fig. 8 of 
Q. While the curve that divides k bins at k = 0.1 Mpc -1 
roughly fits the shape of the systematic feature between 
60 and 180 arcminutes, the curve that uses a division 
at k = 0.01 Mpc -1 is actually a worse fit than the GR 
result. 

However, the main influence leading to our apparent 
detection of a departure from GR is actually due to the 
behavior of the small-angle CFHTLS data (which Q ex- 
cludes out of deference to the uncertainties of non-linear 
modified gravity), as can be seen in Fig. [5] Each of the 
curves in this figure is generated with identical cosmo- 
logical parameters (h = 0.719, Q m = 0.25, and the pri- 
mordial scalar perturbation amplitude, rather than as, is 
also fixed; an exception to this last rule is made for the 
dashed red curve, as discussed in the caption). Post-GR 
parameters are all set to zero except the high k - high z 
value of V — 1 and the high k - low z value of Q — 1, which 
are chosen according to the relationship (which approx- 
imately follows the degeneracy direction of the contours 
drawn in that parameter space) Q—l = — 0.2(V— 1)+0.06. 
One sees that decreasing V — 1 allows the model to bet- 
ter reproduce the precipitous rise of (M^ p ) towards small 
angles. Such small values of V — 1 also reduce the value 
of CTg predicted. The attempt to fit the steep rise in the 
CFHTLS data not only drives as down but this in turn 
then affects other cosmological parameters. 

Figure [6] plots the constraint contours in Q m -as space 
both for post-GR and unmodified GR models. The free- 
dom in the post-GR parameters erases the usual degen- 
eracy between f2 m and as seen in GR, replacing it with 
a degeneracy between as and our post-GR parameters, 
while shifting f2 m . Overall, the MCMC code including 
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FIG. 4: The square of the aperture mass (see Eq. 5 of [i4l |1 
is plotted for different cosmological models in comparison to 
data from the CFHTLS survey. The solid, black curve shows 
the results from the ACDM concordance model in GR. One 
can match the small angle behavior by suppressing growth 
through decreasing the gravitational coupling Q, but increas- 
ing growth by increasing Q m . Exploring the larger angular 
scales, the dashed, red curve shows the effect of changing 
Q(k > 0.01 Mpc -1 ; z < 1) while compensating S7 m . The dot- 
dashed, green curve shows the case for g(k > 0.1 Mpc -1 ; 2 < 
1) as taken in Fig. 8 of [s|. Because this parametrization di- 
vides k bins in the midst of the scales probed by the data, 
this curve fits better the (possibly spurious) bump in (M„ p ) 
seen in the data between 60 arcmin < 6 < 180 arcmin. Data 
is taken from Table B2 of Q. 



CFHTLS data is led to prefer much smaller values of 
a s than are allowed in GR. This extreme shift due to 
the small angle CFHTLS data, and the parametrization- 
dependence exhibited in Fig. [4] due to the CFHTLS data 
bump, give two strong reasons to doubt the significance 
of the exclusion of GR in Fig. [3(d)] Because CFHTLS 
represents the largest current weak lensing data set, we 
continue to use it in the analysis despite these puzzling 
behaviors. However, we will return to these issues in the 
next section and see that COSMOS weak lensing data 
(and CFHTLS data above 10' with regard to the as shift) 
docs not exhibit these deviations. 



III. GALAXY AUTO- AND 
CROSS-CORRELATIONS 

In order to get useful constraints on redshift- and scale- 
dependent deviations from GR, we will need to go beyond 
the basic data sets used so far: WMAP CMB power spec- 
tra 13 1, Union2 supernovae distances 0, and CFHTLS 
weak lensing [l4T |. In particular, different types of cosmo- 
logical probes, more sensitive to density growth, could be 
useful. 

As discussed in Section II of [lj| , the CMB anisotropy 
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FIG. 5: A view of the data in Fig. [4] with a log scale in 6 to 
zoom in on small angles. The theory curves use bins divided 
at fcbin = O.OlMpc" 1 , and each is generated with identical 
background cosmology parameters, fixing the amplitude of 
the primordial scalar perturbations, so that different post-GR 
parameter values give different values of as. Labeled values 
of V — f are set in the high k - high z bin. Values of Q — 1 
in the high k - low z bin are then given by the approximate 
degeneracy relation Q - 1 = -0.2(V - 1) + 0.06. All other 
post-GR parameters are set to zero. To fit the rise at small 
angles, much steeper than in GR, requires very negative V 
and hence low as. Even raising the primordial perturbation 
amplitude (dashed red curve) cannot bring erg into the usual 
range. Values of % 2 reported in the legend are calculated 
naively assuming a diagonal covariance matrix using the error 
bars shown. The four smallest-scale data points are excluded 
from the % 2 calculation. 



spectrum gives poor constraints on modified gravity. 
Note the amorphous, two-lobed shape of the CMB plus 
supcrnovac contours in Figs. 6-9 of that work. This is be- 
cause the ISW term in the CMB auto-correlation goes as 
(4> + ip) and is thus unaware of sign changes induced by 
extreme values of Q and V (or in previous works w). The 
introduction of weak lensing statistics alleviates some of 
this uncertainty. However, much of that ground is lost to 
the introduction of the second post-GR parameter (see 
Fig.®. 

To proceed further, we need to include measurements 
that involve more of the interesting physics of modified 
gravity - further relations between c/>, ip, and A m . The 
two probes we add are the cross-correlation of CMB tem- 
perature fluctuations with the galaxy density field and 
the auto-correlation of the density field, i.e. the galaxy 
power spectrum. 

Section II of [20[ discusses the theory of temperature- 
galaxy cross-correlations. See also Section IV of [2lJ for a 
discussion of how this theory is altered in non-GR grav- 
ity. The salient point is elucidated in Eqs. (4-6) of f20j : 
temperature-galaxy cross-correlations constrain cosmo- 
logical parameters by comparing the matter fluctuations 



FIG. 6: The 68% and 95% cl contours in Q m -as space for 
WMAP7, Union2, and CFHTLS data in the case of our post- 
GR parametrization (solid, green contour) and in the case of 
GR (dashed, blue contour). The inclusion of post-GR param- 
eters seems to eliminate the degeneracy evident in the GR 
case and pulls the contours to lower values of as- This is 
due to the influence of the steeply rising small-scale CFHTLS 
data, as illustrated in Fig. [5] 

A m traced by the galaxy distribution with the sources for 
the metric fluctuations 4> and ip responsible for the ISW 
effect. Because this ISW effect (not its auto-correlation 
in the CMB anisotropy) goes as an integral over redshift 
of 4> + ip times the matter density fluctuation, the cross- 
correlation measurement ends up depending on only one 
factor of cp + ip. Thus, these measurements ought to be 
sensitive to the sign changes that get hidden in the CMB 
anisotropy spectrum. 

It is even likely that temperature-galaxy (Tg) cross- 
correlation data will meaningfully constrain V, since, at 
the small scales considered, the 

V = ^{l + w) 

term in the A m growth equation (A5) of [l[ becomes 
dominant, all other modified gravity terms being sup- 
pressed as (H/k) 2 . 

References [lo, HH provide a module to incorporate 
cross-correlations of the WMAP temperature maps with 
galaxy survey data from the 2-Micron All Sky Survey, 
the Sloan Digital Sky Survey, and the NRAO VLA Sky 
Survey into COSMOMC. We modify this module to ac- 
commodate non-GR values of V and Q and include it into 
our modified COSMOMC. We excise the module code for 
incorporating weak lensing of the CMB described in [23] , 
so as to obtain a clearer picture of the impact of Tg in- 
formation. 

Several works have already applied Tg correlations to 
the ques tion of constraining modified gravity. Refer- 
ence [23| used Tg data to constrain DGP gravity mod- 
els. They found, as in [2(J, that the principal advantage 
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to this data was in constraining models with nonzero 
fifc. Reference [24| considered f(R) gravity and found 
significant improvement over previous constraints using 
just the CMB, though they also found that galaxy clus- 
ter abundances gave constraints that were stronger still. 
These results would seem to indicate that Tg data is not 
as useful at testing gravity as more direct measurements 
of the matter power spectrum. 

However, these studies were carried out in the contexts 
of specific gravity theories in which the relationship be- 
tween high and low k is forced by the theory. Since we 
make no such assumption, we expect (and find) that in- 
clusion of the Tg data significantly improves constraints 
on our high k post-GR parameters. Indeed, [2l[ included 
Tg data in their exploration of a model-independent 
paramctrization based on V (which they call Q) and va 
(their R = 1 + w). Their parameters exhibited a similar 
degeneracy to that discussed in Sec. Ill A[ however the 
linearity of Tg correlations in (0 + VO still allowed them 
to place tighter constraints on the difference V — w than 
CMB and supernova data alone (see their Table 1 and 
Figs. 5-6). 

We also include measurements of the galaxy-galaxy 
(gg) auto-correlation power spectrum of luminous red 
galaxies taken from data release 7 of the Sloan Digital 
Sky Survey and incorporated into COSMOMC by a pub- 
licly available module [25| . These measurements should 
principally be sensitive to V since they are more directly 
measurements of A rn than of <j) + ip, and they are taken 
at scales k > 0.01 Mpc~\ 

At small k, adding temperature-galaxy (Tg) and 
galaxy-galaxy (gg) correlation data produces little 



change in the constraints on V and Q vis-a-vis Figs. 3(a) 
and |3(b)| WMAP constraints from the ISW effect domi- 
nate at low k, plus there are no gala xy-g alaxy data points 
at k < 0.03 h Mpc -1 (see Fig. 8 of [$5|). 

For large k 7 however, where the CMB provides gener- 
ally poor constraints, the addition of temperature-galaxy 
and galaxy-galaxy data can significantly alter limits on 
our post-GR parameters. Figure [7] illustrates two ex- 



amples of this. Figure 7(a) shows the strengthening of 



constraints in the 95% cl contour from Fig. 3(c) upon 
incorporating as well the temperature-galaxy data, and 
the Tg plus galaxy-galaxy data. With current galaxy- 
galaxy data, most of the improvement is due to Tg, but 
one can anticipate that as larger galaxy surveys includ- 
ing next generation surveys are completed then galaxy 
power spectra will become an important ingredient in 
testing gravity (see Sec. |IV] for future projections). In 
particular, V is still not well determined now. 

Conversely, Fig. |7(b)| shows the effects of adding Tg, 
and then Tg plus gg, can for some variables shift the 
contours instead of tightening them. This may represent 
a certain tension between data sets; it is interesting to 
note that Fig. |7(b)| finds deviation from GR at the 95% 
cl, and we return to the role of CFHTLS tension in this 
below. 

Figures [5] update all of the plots in Figs. [3] using all of 
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FIG. 7: [Top panel] 95% cl contours in V-Q space for the high 
k, high z bin are compared for three different combinations 
of data sets. The solid, black contour shows the results from 
Fig. |3(c)] using CMB, supernovae, and weak lensing data. The 
thin-dashed, red contour adds temperature-galaxy (Tg) cross- 
correlation data from [2(|. The thick-dashed, green contour 
further adds galaxy-galaxy (gg) correlation data from [25l [. 
The diagonal, dot-dashed line gives the n = boundary. As 
datasets are added, the contours close in on GR parameter 
values (the magenta x). Current galaxy correlation data is not 
yet sensitive enough though to put a meaningful constraint on 
V. [Bottom panel] Addition of data sets can sometimes shift 
rather than tighten the contours, as shown here in Q(l < 
z < 2)-Q(z < 1) space for the low k bin. Note that with the 
additional data GR now lies on the edge of the 95% cl region. 



the data sets discussed. The results are the foreground, 
blue contours. We see that Q is constrained with an un- 
certainty of roughly 0.1, while V is unknown to within 
~ 1. All cases, except high k - high z, have pulled further 
off the fj, — restricted area. Note that the 95% cl con- 
tour in the high k - low z bin, Fig. |8(d)| using CFHTLS 
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FIG. 8: 68% and 95% cl constraints on V — 1 and Q — 1 are plotted for the two redshift and two wavenumber bins, using CMB, 
supernovae, weak lensing, Tg, and gg data. Foreground, blue contours use CFHTLS weak lensing data. Background, yellow 
contours use COSMOS weak lensing data. The dotted contours reproduce the 95% cl contours without Tg or gg data from 
Figs. [21 The diagonal, dot-dashed line gives the /j, = boundary, from which the low k contours at least have now pulled away. 
The x's denote GR values. Both k bins at low z exhibit some preference for non-GR parameter values when using CFHTLS, 
but not when using COSMOS, weak lensing data. 



data, still excludes General Relativity, although as we 
have stated this is possibly due to systematics in the 
CFHTLS weak lensing data. It now appears that the low 
k - low z bin also prefers non-GR values of our parame- 
ters, though in this case the apparent exclusion of GR is 
just at 95% cl. Since this effect did not manifest itself un- 
til we added the galaxy-based datasets, this could either 
be an effect of systematic tension between galaxy-count 
measurements and other data sets, or a true restriction 
from the increased precision. Note that Fig. |7(b)| gives 
another view of the low k - low z deviation in Q. 

To test the hypothesis that the exclusion of GR at low 
z is due to systematic effects in the CFHTLS data, we 
plot the same constraints substituting weak lensing data 



from the COSMOS survey in the place of CFHTLS 
data. COSMOS constraints are the background, yellow 
contours in Fig. [5J The constraints in the low k bins ap- 
pear almost unaffected by this substitution (though the 
low k - low z bin no longer hints at an exclusion of GR, 
as it did in the case of the CFHTLS data) . This should 
not be surprising. The low k bins correspond to scales 
where the WMAP7 data has a lot of constraining power. 
The high k high z bin also appears moderately insen- 
sitive to which weak lensing set is used. In the high k 
- low z bin, though, we find that the 99% cl exclusion 
of GR vanishes when COSMOS is used, and GR instead 
lies comfortably within the 68% cl contour. This could 
mean that the COSMOS data is less subject to spuri- 
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FIG. 9: 68% and 95% cl constraints on the correlations of V — 
1 (top panel) and Q — 1 (bottom panel) between redshift bins, 
using CMB, supernovae, weak lensing (CFHTLS), Tg, and gg 
data. Contours are labeled according to k binning. The x's 
denote GR values. The high k case of V exhibits a deviation 
from GR, corresponding to a different growth amplitude. 



ous systematic effects. Note that the combination using 
COSMOS data is somewhat less constraining due to the 
small sky area of COSMOS. It will be interesting to see 
what results occur once we have data from larger, more 
detailed future weak lensing surveys. 

We can also reexamine the redshift dependence of each 
of the post-GR parameters. Figures [9] plot the correla- 
tions between redshift bins for V and Q, for each bin of 
wavenumber k using CFHTLS data. The values of Q fol- 
low k are positively correlated, as found before for scale- 
independent Q, while for high k the values are rather in- 
dependent. This can be understood by considering that 
the small k bin has the greatest effect on the ISW im- 
print in the CMB. Since the ISW effect goes as (<fi + ip) 2 , 



the CMB data will prefer parameter combinations that 
minimize the change in Q across redshift bins. Because 
the large k bin encompasses scales over which the ISW 
effect is subdominant, such a preference is not operative 
there. 

For V there is negative correlation as before, for both 
low and high k. This is a manifestation of the role of V 
in regulating the growth of structure, i.e. A m . The low 
k contours in Fig. 



have a main degeneracy direction 
parallel to lines with a slope in {V(l < z < 2), V(z < 1)} 
space of ~ —1.3. The high k contours have a degener- 
acy direction of approximate slope <~ —1.2. Integrating 
Eq. (A5) of Q 

for values of O m close to those favored by 
WMAP and Union2 (Q m ~ 0.25), one finds that (inde- 
pendent of Q and H ) values of {V(l < z < 2), V(z < 1)} 
that lie along lines with slopes ranging between —1.3 
and —1.4 return values of the relative growth A m (z = 
0)/A m (z = 100) that are similar to within a few percent. 
Displacement perpendicular to this direction controls the 
absolute growth factor. The offset of the high k contours 
from the low k contours (and from GR), signifies a pref- 
erence for suppressed growth relative to GR in the high 
k modes; GR lies outside the 95% cl contour in the high 
k bin. Again, this could be due to the odd bump in the 
CFHTLS weak lensing power or the steep rise towards 
small angles in Fig. [5] Figures [TU] show the effect of us- 
ing COSMOS weak lensing data instead of CFHTLS data 
on the correlations of V — 1 and Q — l across redshift bins. 
As before, the exclusion of GR vanishes. 

As in the scale-independent case, we find that our post- 
GR parameters correlate most strongly with as out of all 
of the usual cosmological parameters. Once again, larger 
V amplifies growth and induces a larger a$ while larger 
Q brings lower values of <j% into agreement with the data. 
These correlations are only manifest in the high k bin, in- 
dicating that they are principally dependent on the weak 
lensing, Tg, and gg data sets. The correlation with V ap- 
pears in the high z bin, as expected due to the cumulative 
effect of growth over time (see discussion in Sec. Ill C[) . In 
the case of Q, probing the potentials, the correlation with 
as appears in the low z bin, as expected for the weak 
lensing data weighted toward z < 1. The shift to low as 
as seen in Fig. [5] does not occur when COSMOS data is 
used, or when CFHTLS data is restricted to 9 > 10'. 



IV. CONSTRAINTS POSSIBLE WITH FUTURE 
DATA 

Given the weak constraints on V in Fig. [SI and possible 
hints of deviations from GR, it is important to investigate 
the capabilities of future galaxy surveys. These should 
provide us with more direct measurements of, and much 
better precision on, the growth of density perturbations 
through the galaxy power spectrum, testing GR and im- 
proving our knowledge of post-GR parameters. 

We consider the specific example of BigBOSS [26jj, 
a proposed ground-based survey intended to constrain 
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cosmology by measuring the baryon acoustic oscillations 
and redshi ft sp ace distortions in the galaxy distribution. 
Reference [2Zj explored BigBOSS tests of gravity (and 
dark sect or p hysics) in terms of the gravitational growth 
index 7 [10|, using a Fisher matrix calculation. Here, 
we carry out a more sophisticated Markov Chain Monte 
Carlo fit to simulated data, allowing for scale- and time- 
dependence in the gravitational modifications through 
our binned k, binned z post-GR parameters. 

We generate our mock BigBOSS data around the 
ACDM, GR (V = Q = 1) maximum likelihood cosmology 
of WMAP7. The data is considered in the form of mea- 
surements of the redshift-distorted galaxy-galaxy power 
spectrum 



1800 supernovae out to z — 1.5 ("JDEM") in our pro- 
jected data MCMC calculation. For computational effi- 
ciency we do not consider the gravitational lensing of the 
CMB. 

FigureslTTIshow the 68% and 95% confidence limit con- 
tours resulting from the full MCMC calculation on our 
mock data sets. There is an order of magnitude or more 
improvement in the constrai nt p laced on V in all bins by 
including BigBOSS. As in df, we find that BigBOSS 
will be able to constrain departures from the growth his- 
tory of GR (here parametrized as V — 1, there by 7) to 
within ~ 10%. Constraints on Q also improve, though by 
more modest factors. 



P g (k; z) = P g {k, n; z) = (b+ ffi 2 ) 2 P(k; z) 



(11) 



where P is the power spectrum of matter ovcrdensities, b 
is the linear bias relating the overdensity of galaxy counts 
to the overdensity of matter, / = <iln A m (fc, a) /din a is 
the growth factor, and /1 is the cosine of the angle k makes 
with the line of sight. 

We take surv ey p arameters, including galaxy number 
densities, from [28|, and consider emission line galaxies 
(ELG) and luminous red galaxies (LRG) as two separate 
data sets. The bias & is a function of redshift, 



b(z) = b Q 



A ?ra (fc,z = 0) 
A TO (fc, z) 



(12) 



where &o is a nuisance parameter to be marginalized over 
for each data set. The fiducial values are &o,elg = 0.8 
and 60. lrg = 1-7. The values and redshift dependence 
are good fits to current galaxy observations, and can be 
motivated by comoving clustering models [29J. 

In calculating the galaxy power spectrum we consider 
modes 10~ 4 < k < 10 -1 Mpc _1 . As in [30j], we assume 
that the covariance matrix of P g is diagonal with 



1C 3r r Pg(ki, /J*m)Pg(kj, Mn) 
1D7T O l jO mn 



1 + nP g (ki,n 7 



V d 3 k 
\ i 2 



(13) 



where Vb is the real space volume of the survey and n is 
the selection function of the survey as a function of co- 
moving distance r. Thus, the likelihood C for a cosmolog- 
ical model which predicts galaxy-galaxy power spectrum 
P g {k,n) is given by 



-21n[£] = 



1 

4tt 



A dr 2 

dz—r 
dz 



d\nkk J 



(Pg - Pg? ( n(r)P g 



pi 
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(14) 



1 + n(r)P g , 

We also include mock Planck CMB dat a g enerated 



with the COSMOMC module FuturCMB 31, \M and 



mock future supernova data based on a space survey of 



V. CONCLUSIONS 

The suite of current cosmological data has grown to the 
point that increasingly sophisticated model-independent 
extensions to general relativity can be tested. This in- 
cludes both time- and scale-dependent modifications; we 
utilize bins in redshift z and wavemode k for localization 
of the effects and clarity of physical interpretation. The 
functions Q and V investigated here in detail, giving a 
complete model-independent description (together with 
stress-energy conservation) of the gravitational modifica- 
tions, are closely tied to the sum of the metric potentials 
and to the matter growth, respectively. They also have 
the virtue of being substantially decorrelated from each 
other. On the other hand, correlations across redshift or 
across length scales can be easily studied, giving deeper 
insight into the effects of the modifications and where 
they show up in the observations. 

Using current CMB, supernova, and weak lensing data 
from CFHTLS we find an inconsistency with general rela- 
tivity at near the 99% confidence limit at k > 0.01 Mpc -1 
and z < 1 . Through a series of investigations we identify 
its origin as being due to an abnormally steep rise in the 
weak lensing power at small scales. This rise strongly 
shifts V from the GR value and also drives down the 
estimated value of cr$. The CFHTLS data also shows 
an unusual bump in the power at larger scales. By re- 
placing the CFHTLS measurements with COSMOS weak 
lensing data, we find that all these deviations vanish and 
that GR provides an excellent fit. Thus, the deviations 
may originate in systematic effects in the CFHTLS data 
(a possibility also raised by members of the CFHTLS 
team) interacting with increased freedom from the post- 
GR parameter fitting. 

The addition of galaxy clustering measurements, 
through both the CMB temperature-galaxy count (Tg) 
and galaxy-galaxy power spectrum (gg) statistics, tight- 
ens the constraints on smaller scales (high k). This im- 
provement is especially noticeable in V, since it enters 
in the matter growth. Again, the full combination us- 
ing CFHTLS data shows significant deviations, which go 
away on substitution with COSMOS data. 

Given the important role of galaxy survey data, and 
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the still weak constraints on the V deviation parameter 
(only of order unity), we examine the potential leverage 
of future galaxy survey measurements, specifically from 
BigBOSS. These appear quite promising for confronting 
general relativity with further measurements, giving a di- 
rect probe of growth and one that could be highly precise 
from the large statistics. Together with Planck CMB and 
supernovac data, such a galaxy survey could improve the 
area uncertainty on the post-GR parameters in each of 
four bins of redshift-scale by factors from 10 to 100. This 
is an exciting prospect as we seek to understand gravity 
as the most pervasive and dominant force in the universe. 
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FIG. 10: 68% and 95% cl constraints on the correlations 
of V — 1 and Q — 1 between redshift bins, using CMB, su- 
pernovae, weak lensing, Tg, and gg data. Substituting COS- 
MOS (background) for CFHTLS (foreground) eliminates the 
apparent exclusion of GR. The x's denote GR values. We do 
not plot the high k bin correlations for Q — 1 as there are no 
interesting correlations or deviations there. 
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FIG. 11: 68% and 95% cl constraints on V — 1 and Q — 1 are plotted for the two redshift and two wavenumber bins using 
mock future data. Foreground (blue) contours use mock BigBOSS, Planck, and JDEM supernova data. Background (yellow) 
contours use only mock Planck and JDEM supernova data. The dotted contours recreate the 95% cl current data contours 
from Figs. [H] (using CFHTLS) to illustrate the expected improvement in constraints. The x's denote the fiducial GR values. 



